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Abstract. This paper is concerned with the optimization of 
the trajectories of an aerospace plane i This is a hypervelocity 
vehicle capable of achieving orbital speed, while taking off 
horizontally. The vehicle is propelled by four types of engines: 
turbojet engines for flight at subsonic speeds/ low supersonic 
speeds; ramjet engines for flight at moderate supersonic 
speeds/low hypersonic speeds; scramjet engines for flight at 
hypersonic speeds; and rocket engines for flight at near— orbital 

speeds . 

A s ingle- stage-to-orbit (SSTO) configuration is considered, 
and the transition from low supersonic speeds to orbital speeds 
is studied under the following assumptions: the turbojet portion 
of the trajectory has been completed; the aerospace plane is 
controlled via the angle of attack a(t) and the power setting 
3 (t) ; the aerodynamic model is the generic hypersonic aerodynamics 
model example (GHAME) . Concerning the engine model, three options 
are considered: (EMI) this is a ramjet/scramjet combination in 

which the scramjet specific impulse tends to a nearly-constant 
value at large Mach numbers; (EM2) this is a ramjet/scramjet 
combination in which the scramjet specific impulse decreases 
monoton ic al ly 1 at large Mach numbers; (EM3) this is a 
ram jet/ scram jet/rocket combination in which, owing to stagnation 
temperature limitations, the scramjet operates only at M <_ 15 ; at 
higher Mach numbers, the scramjet is shut off and the aerospace 
plane is driven only by the rocket engines. 
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Under the above assumptions, four optimization problems are 
solved using the sequential gradient-restoration algorithm for 
optimal control problems: (PI) minimization of the weight of fuel 

consumed; (P2) minimization of the peak dynamic pressure; (P3) 
minimization of the peak heating rate; and (P4) minimization of 
the peak tangential acceleration. 1 The above optimization studies 
are carried out for different combinations of constraints, 
specifically: initial path inclination either free or given 
(Yq = o) ; dynamic pressure either free or bounded (q £ 1500 lbf/ft ) ; 
tangential acceleration either free or bounded (a T <_ 3g e ) . 

The preliminary conclusions are as follows : 

(a) For an aerospace plane governed by GHAME + EMI, the 
SSTO mission requires a weight of fuel consumed equal to 34.3% of 
the initial weight. 

(b) For an aerospace plane governed by GHAME + EM2 , the 
SSTO mission requires a weight of fuel consumed equal to 44.3% of 
the initial weight. 

(c) For an aerospace plane governed by GHAME + EM3 , the 
SSTO mission requires a weight of fuel consumed equal to 60.7% of 
the initial weight. 

(d) If one assumes that engine model EM2 is the one closer 
to reality, then the SSTO mission appears to be feasible. 

Obviously, its ability to deliver payloads can be improved via 
progress in the areas of aerodynamic properties and specific 
impulse properties. 
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(e) If one assumes that engine model EM3 is the one closer 
to reality, then the SSTO mission appears to be marginal, unless 
substantial progress is achieved in the areas of aerodynamic 
properties and specific impulse properties. Under this scenario, 
alternative consideration should be given to studying the 
feasibility of a two-stage-to-orbit (TSTO) mission. 

Key Words . Flight mechanics, hypervelocity flight, 
atmospheric flight, optimal trajectories, aerospace plane, 
sequential gradient-restoration algorithm. 
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1. Introduction 

The aerospace plane is a hypervelocity aircraft which must 
take off horizontally, achieve orbital speed, and then land 
horizontally. At this time, its configuration is not precisely 
known, but it can be assumed that the powerplant includes the 
combination of four types of engines: turbojet engines for flight 
at subsonic speeds/ low supersonic speeds; ramjet engines for 
flight at moderate supersonic speeds/ low hypersonic speeds? 
scramjet engines for flight at hypersonic speeds; and rocket 
engines for flight at near-orbital speeds. 

In this paper, we refer to a single-stage-to-orbit (SSTO) 
configuration and we study the transition from low supersonic 
speeds to orbital speeds under the following assumptions: (i) the 

turbojet portion of the trajectory has been completed; (ii) the 
aerospace plane is controlled via the angle of attack ot(t) and 
the power setting $(t) ; (iii) the switch times from one powerplant 
to another are parameters being optimized. 

Concerning the aerodynamics model, two configurations have 
been considered thus far in the aerospace plane literature: the 
generic hypersonic aerodynamics model example (GHAME) and the 
Langley accelerator model example (LAME) . The first of these 
configurations (GHAME) is considered here. 

Concerning the engine model, three options are considered: 
(EMI) this is a ramjet/scramjet combination in which the 
scramjet specific impulse tends to a nearly-constant value at 
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large Mach numbers; (EM2) this is a ramjet/scramjet combination 
in which the scramjet specific impulse decreases monotonically at 
large Mach numbers; (EM3) this is a ramjet/scramjet/rocket combination 
in which, owing to stagnation temperature limitations, the scramjet 
operates only at M <_ M*; at higher Mach numbers, the scramjet is 
shut off and the aerospace plane is driven only by the rocket engines. 
Here, M* is a threshold Mach number. 

With the above understanding, we study four basic 
optimization problems: (Pi) minimization of the weight of fuel 

consumed; (P2) minimization of the peak dynamic pressure; (P3) 
minimization of the peak heating rate; and (P4) minimization of 
the peak tangential acceleration. These optimization problems are 
solved by means of the sequential gradient-restoration algorithm (SGRA) 
for different combinations of constraints imposed on the initial 
path inclination y q , the dynamic pressure q, and the tangential 
acceleration a T . Specifically ,y 0 can either be free or given 
(Yq = 0); q can either be free or bounded (q £ 1500 lbf/ft^);and a T 
can either be free or bounded (a T £ 3g e ) . 

Previous studies of interest for the aerospace plane can be 
found in Refs. 1-19. References 1-3 contain basic concepts. 

References 4—8 deal with flight mechanics, hypervelocity flight, 
and propulsion systems. References 9—12 treat design problems, 
while Refs. 13-18 consider trajectory optimization and guidance. 
Concerning aerodynamics, two widely used configurations are the 
generic hypersonic aerodynamics model example (GHAME, Ref. 19) 
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and the Langley accelerator model example (LAME) . While GHAME has 
been used in Ref. 15, LAME has been considered in Ref. 16. 
Concerning the engine models, we note that EMI has been 
considered in Ref. 15, while EM2 has been considered in Refs. 
16-18. A model akin to EM 3 has been considered in Ref. 10, albeit 
with a lower threshold Mach number than that considered in this 
paper . 

Section 2 contains the notations, and Section 3 presents the 
system description. Sections 4—5 deal with the performance 
indexes and present a classification of the problems being 
studied. In Section 6, we present the results obtained on optimal 
trajectories for engine models EMI, EM2 , EM3 . Finally, the 
conclusions are given in Section 7 . 
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2 . Notations 

Throughout this paper, the following notations are employed: 
a = acceleration, ft/sec 2 ; 

C D = drag coefficient; 

C L = lift coefficient; 

D = drag, lbf; 

E - lift-to-drag ratio; 

g = local acceleration of gravity, ft/sec^; 
g e = sea-level acceleration of gravity, f t/ sec z ; 
h = altitude, ft; 

I S p - specific impulse, sec; 

L = lift, lbf; 
m = mass, lbf sec 2 /ft; 

M = Mach number; 
q = dynamic pressure, lbf/ ft 2 ; 

Q = heating rate, BTU/ft 2 sec; 

r = radial distance from the center of the Earth, ft; 

r e = radius of the Earth, ft; 

S = reference surface area, ft^; 

p 

S e = combustor cross-sectional area, ft^; 
t = dimensionless time; 

T = thrust, lbf; 

V = velocity, ft/sec; 

W = mg e = sea-level weight, lbf; 
x = distance along the Earth surface, ft; 

= angle of attack, rad; 


a 



5 


AAR-247 


3 = power setting; 

Y = path inclination, rad; 

6 as inclination of the thrust with respect 
to the aircraft reference line, rad; 

6 = running time, sec; 

y = Earth's gravitational constant, ft 3 /sec 2 ; 
p = air density, lbf sec 2 /ft 4 ; 

t = final time, sec. 

Subscripts (EMI + EM2) 

0 = beginning of ramjet phase/initial point; 

1 = beginning of scram jet phase; 

2 = end of scramjet phase/ final point. 

Subscripts (EM3) 

0 = beginning of ramjet phase/initial point; 

1 = beginning of scramjet phase; 

2 = beginning of rocket phase; 

3 = end of rocket phase/final point. 

Superscript 

• = derivative with respect to dimensionless time. 

Acronyms 

GHAME ~ general hypersonic aerodynamics model example; 
LAME = Langley accelerator model example; 

SGRA = sequential gradient-restoration algorithm; 

SSTO = single-stage-to-orbit; 

TSTO = two-stage-to-orbit . 
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3 . System Description 

We consider a single-stage-to-orbit (SSTO) aerospace plane, 
powered by the combination of turbojet/ramjet/scramjet/rocket 
engines. We employ the following hypotheses: (i) the turbojet 

portion of the trajectory has been completed; (ii) flight takes 
place in a vertical plane over a spherical Earth; (iii) the Earth's 
rotation is neglected; (iv) the gravitational field is central 
and obeys the inverse square law; (v) the aerospace plane is 
controlled via the angle of attack a(t) and the power setting 3 ( t ) - 

3.1. Time Normalization . We denote with 0 the actual time 
and with t the normalized time. The normalization is done in such 
a way that the normalized time duration of each segment of the 
trajectory is one. Hence, the transformation relations are as follows 


6 = T 1 t, 

0 

< t 

< 

1 

(ramjet) , 

(la) 

0 = T 1 + X 2 (t-1) , 

1 

1 A 
f+ 

< 

2 

(scram jet) , 

(lb) 

0 = T, + Try + T_(t“2), 

2 

< t 

< 

3 

(rocket) , 

(lc) 


with the implication that 


0 = T x , 

0 

< 

t <_ 

1 

(ramjet) , 

(2a) 

* 

0 = t 2 . 

1 

< 

t < 

2 

(scramjet) , 

(2b) 

• 

0 = T-, 

2 

< 

t < 

3 

(rocket) , 

(2c) 


and that 
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0 O = °' 


(3a) 


9 1 = V 


(3b) 


9 2 T 1 + T 2 ' 


(3c) 


0 3 T 1 + T 2 + t 3 * 


(3d) 


Note that is the final time for engine models EMI, EM2 and 
that 0^ is the final time for engine model EM3 . 

3.2. Differential System. With the above assumptions and 
upon normalizing the time duration of each segment of the trajectory 
to unity, the motion of the aerospace plane is described by the 
following differential system: 


x = t [ ( r g /r) Vcosy] , (4a) 

h = t [V siny] , { 4b) 

V = t [ (Tg /W)cos(a + 6) - Dg /W - gsiny] , (4c) 

e e 

y = x t (Tg^/WV) sin (a +6) + Lg e /WV + (V/r - g/V)cosy], (4d) 

W = t [ “T/I ] . (4e) 

sp 


Here, the dot denotes derivative with respect to the normalized 
time and the flight duration t takes the following values: T = T ]_ 
for the ramjet segment, t — T 2 for the scramjet segment, and 
t — t ^ for the rocket segment of the trajectory. 

In the above system, the following functional relations hold: 


= r + h, 
e 


r 


(5a) 
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g = y/r 2 = y/(r + h) 2 , (5b 

2 

where y = g e ^ e denotes the Earth's gravitational constant. 

The weight W = mg e appearing in Eqs. (4) is the so-called sea-level 
weight, which is based on the sea-level acceleration of gravity 
g . Such weight differs from the local weight mg, which is based 
on the local acceleration of gravity g. 

3.3. Aerodynamic Data. The drag and the lift are given by 


D = <l/2)C D pSV , 


L = <l/2)C L pSV , 

with p = p (h) . Generally speaking, the aerodynamic coefficients 
C D , C L depend on the angle of attack a, the Mach number M, and 
the Reynolds number Re. If the dependence on the Reynolds number 
is disregarded, the aerodynamic coefficients take the form 


C D = C D (a,M) , 


C L = C L (a,M) , 

.th the implication that the lift-to-drag ratio E = L/D - C L /C D 


is a function of the form 


E = E (a , M) . 


The functions (7) are plotted in Fig. 1 with reference to 
the generic hypersonic aerodynamics model example (GHAME) . 

For computational purposes, it is convenient to approximate 
the aerodynamic coefficients with polynomial relations of the type 
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C D - A q (M) + A 1 (M)a + A 2 (M)a 2 , (8a) 

2 

C L = B q (M) + B 1 (M)a + B 2 (M)ct , (8b) 

with the implication that 

2 2 

E = [B q (M) + B x (M)a + B 2 (M)a ]/[A Q (M) + A x (M) a + A 2 (M) a ]. (8c) 

The coefficients A^ (M) , B^ (M) are computed by means of a least 
square fit of the available GHAME data at various Mach numbers 
and angles of attack. 

3.4. Engine Data . For the ramjet engines, the following 
simplified representation is assumed for the thrust and the specific 
impulse : 


T = (M) p/p* , 


( 9a) 


I 

sp 


I (M) , 
sp* x ' ’ 


(9b) 


with the implication that the fuel rate (weight of fuel consumed 
per unit time) is given by 


T/I = $[T ( M ) / I (M) ] p/p* . 

sp * sp* ' 


(9c) 


Here, p* is a reference density (density at the reference altitude 
h*= 100 kft) , T* (M) is a reference thrust (thrust for 3-1 and 
h = h*) , and I g p*(M) is a reference specific impulse (specific 
impulse for 3=1 and h = h*) . While the thrust is assumed to 
depend on the power setting, the altitude, and the Mach number, 
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the specific impulse is assumed to depend only on the Mach number; 
the dependence of the specific impulse on the power setting and 
the altitude is disregarded in line with the feasibility character 
of the present study. For the same reason, the dependence of the 
thrust and the specific impulse on the angle of attack, relevant 
to a precision study, is disregarded within the bounds of the 
present feasibility study. 

For the scramjet engines, the representation (9) is retained. 

However, the reference thrust T* (M) and the reference specific 

impulse I (M) are now described by different functions, 
sp* 

For the rocket engines, the following simplified relations 
are assumed: 


T = 3T*, 


(10a) 


I 


sp 


sp* ' 


(10b) 


with the implication that the fuel rate is given by 


T/I sp = 3(T*/I spi ). (10c) 

Here, T. is a reference thrust (thrust for 3=1) and I is 

sp* 

a reference specific impulse (specific impulse for 6=1). Both 

T. and I are assumed to be constant. This means that the weak 
* sp* 

dependence of these quantities on the altitude is disregarded. 

In this paper, three engine models are considered: 


(EMI) This is a ram j et/scram j et combination in which the 
scramjet specific impulse tends to a nearly-constant value at 
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large Mach numbers; see Ref. 15. For this combination, the functions 
(9) are shown in Fig. 2 (ramjet) and Fig. 3 (scramjet) under the 

p 

assumption that the combustor cross-sectional area is S = 400 ft . 

e 

(EM2) This is a ram j et/scramj et combination in which the 
scramjet specific impulse decreases monotonically at large Mach 
numbers; see Refs. 16-18. For this combination, the functions (9) 
are shown in Fig. 2 (ramjet) and Fig. 4 (scramjet) under the 

2 

assumption that the combustor cross-section area is S = 400 ft . 

e 

(EM3) This is a ramj et/scramjet/rocket combination in which, 

owing to stagnation temperature limitations, the scramjet operates 

only at M <_ 15; at higher Mach numbers, the scramjet is shut off 

and the aerospace plane is driven only by the rocket engines.. For 

this combination, the functions (9) are shown in Fig. 2 (ramjet) 

and Fig. 5 (scramjet, M <_ 15) under the assumption that the 

2 

combustor cross-sectional area is S e = 400 ft ; the functions (10) 

are shown in Fig. 6 (rocket, M > 15) under the assumption that 

the reference thrust (hence, the maximum thrust) is T* = 189200 Ibf . 

3.5. Control Inequality Constraints . To obtain realistic 
solutions, the presence of upper and lower bounds on the angle 
of attack and the power setting is necessary. Therefore, the 
two-sided inequality constraints 

a n < a < a , (11a) 

% — — u 

h < ^ < e u (llb) 

must be satisfied everywhere along the interval of integration. 
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The inequality constraints (11) can be converted into equality 
constraints by means of trigonometric transformations of the 
type 


a = (1/2) (a • + a ) + (1/2) (a^ - a £ )sinu, (12a) 

3 = (1/2) (3^ + 3 U ) + (1/2) (3 u - 3^)sinw. (12b) 

Therefore, the angle of attack a(t) is replaced with the auxiliary 
control u(t), while the power setting 3(t) is replaced with the 
auxiliary control w(t) . After a solution is found for the auxiliary 
controls, then the original controls are computed with (12). 

For the GHAME configuration, the bounds (11a) are given by 

- -2.0 deg, = 12.0 deg. (13a) 

These bounds are dictated solely by the availability of data. For the 
engine models EMI, EM2 , EM3 , the bounds (lib) are given by 

- 0, 3 U = 1. (13b) 

3.6. Derived Quantities . After a solution of Eqs . (4) 

is available, certain derived quantities can be computed. The more 
relevant ones are listed below: 

(i) the dynamic pressure 

q = (1/2) pV 2 ; (14a) 

(ii) the heating rate 


Q = C / ( p/p * ) (V/V*) 


3.07 


(14b) 
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here, p* is a reference density (density at h*= 100 kft) and V* 

is a reference velocity (V* = 10 kft/sec) ; under the assumption that 

the nose radius is r = 1.0 ft, the constant C has the value 

n 

2 

C= 101.92 BTU/ft sec; 

(iii) the tangential acceleration 


a T - (Tg e /W)cos(a +6) - Dg e /W - gsiny; 


(14c) 


(iv) the normal acceleration 

2 

a N - (Tg e /W)sin(a + 6) + Lg e /W + (V /r - g)cosy; (14d) 

(v) the total acceleration 

a = /(a 2 + a 2 ) . (14e) 

3.7. Supplementary Bounds . In addition to the control 
inequality constraints (11), supplementary bounds can be imposed 
on the quantities (14). For instance, a bound of the form [see (14a)] 

q <: q u - 1500 lbf/ft 2 (15a) 

is a state constraint, while a bound of the form [see (14c)] 


T 


< a 


Tu 


= 3g 


(15b) 


is a state/control constraint. 

The simplest way to account for the bounds (15) is by means 
of penalization techniques. The functional being minimized 
(see Section 4) is augmented by a penalization functional with 
integrand 
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C l (q ‘ q u )n + C 2 (a T " a Tu )R ' (15C) 

where C^, C 2 are constants and n is a suitable integer. The 
value of the first constant is ^ = 0 if (15a) is satisfied and 
C 1 > 0 if (15a) is violated; the value of the second constant 
is = 0 if (15b) is satisfied and C 2 > 0 if (15b) is violated; 
the desirable value of the integer is n = 3, since this ensures 
the continuity of the first and second derivatives of the integrand 
of the penalization functional at points located on the constraint 
boundaries . 

3.8. Boundary Conditions . For engine models EMI, EM2 , EM3 , 
the dimensionless time t = 0 marks the end of the turbojet phase, 
the beginning of the ramjet phase, as well as the initial time. 

At t = 0 , we assume that 


x o = 0 ft ' 

h Q = 42004 ft = 12.8 km, 

V Q = 1936 ft/sec, 

Yq = free or = 0.0 deg, 
W Q = 290000 lbf , 
with the implication that 
M q = 2, 

q 0 = 1000 lbf /ft 2 . 


(16a) 

(16b) 

(16c) 

(16d) 

(16e) 

(16f) 


(16g) 
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For engine models EMI, EM2 , EM3 , the dimensionless time 
t = 1 marks the end of the ramjet phase and the beginning of the 
scramjet phase. At t = 1, the continuity of all the state variables 
is required. 

For engine models EMI and EM2 , the dimensionless time t = 2 
marks the end of the scramjet phase as well as final time. At 


t = 2, we assume that 

x 2 = free, (17a) 

h 2 = 262467 ft = 80.0 km, (17b) 

V 2 = 25792 ft/sec, (17c) 

Y 2 = 0.0 deg, (17d) 

W 2 = free, (17e) 

with the implication that 

M 2 = 27.8, ( 17f ) 

q 2 = 11.9 lbf/ft 2 , ( 17g) 


and that orbital speed is achieved. 

For engine model EM3 , the dimensionless time t = 2 marks the 
end of the scramjet phase and the beginning of the rocket phase. 

At t = 2, we assume that 

M 2 = 15. (18) 

In addition, we require the continuity of all the state variables. 
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For engine model EM3 , the dimensionless time t = 3 marks 
the end of the rocket phase as well as the final time. At t = 3 , 


we assume that 

= free, (19a) 

h 3 = 262467 ft = 80.0 km, (19b) 

V 3 = 25792 ft/sec, (19c) 

Y 3 = 0.0 deg, (19d) 

W 3 = free, (19e) 

with the implication that 

M 3 = 27.8, (19f ) 

q 3 = 11.9 lbf/ft 2 , (19g) 


and that orbital speed has been achieved. 

3.9. Summary . The relations governing the motion of the 
aerospace plane include: the differential system (4); the control 
inequality constraints (11) , converted into control equality 
constraints by means of the trigonometric transformations (12) ;the 
possible presence of the supplementary bounds (15) , accounted via 
penalization techniques; the boundary conditions (16) -(17) for 
engine models EMl , EM2 ; and the boundary conditions (16), (18), 

(19) for engine model EM3 . 

In this formulation, the independent variable is the time t, 
which varies in the range 0 <_ t <_ 2 for engine models EMI, EM2 
and in the range 0 < t < 3 for engine model EM3 . The dependent 
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variables include five state variables [x(t), h(t), V(t), y(t), W(t)], 

two control variables [u(t), w(t)], plus two parameters [x^, t 2 ] 

for engine models EMI, EM 2 and three parameters [t^, t 2 , t^] for 

engine model EM3 . After a solution is found for the auxiliary 

control variables, the original control variables [a(t), B(t)3 

are recovered via the trigonometric transformations (12) . 

3.10. Experimental Data . The following data are used in the 

numerical experiments on optimal trajectories. 

Spaceplane . For the aerospace plane, the initial weight 

(weight at the end of the turbojet phase) is Wq = 290000 lbf; 

2 

the reference surface area (wing area) is S = 6000 ft ; the 

aerodynamic data for the GHAME configuration are given in Fig. 1; 

the angle of attack is subject to the inequality -2.0 <_ a ^12.0 deg. 

Engines. The data for engine models EMI, EM2 , EM3 are given 

in Figs. 2-6. For all models, the combustor cross-sectional area 

2 

of both the ramjet and the scramjet is S e = 400 ft . For engine 

model EM3 , the maximum rocket thrust is T* = 189200 lbf and the 

rocket specific impulse is I = 444 sec. The inclination of the 

thrust with respect to the aircraft reference line is 6 = 0.0 deg; 

the power setting is subject to the inequality 0 £ 3 £ 1. 

Physical Constants . The radius of the Earth is assumed to be 

r = 0.2093E+08 ft = 6378 km. The Earth’s gravitational constant 
e 

3 2 . 

is y = 0.1409E+17 ft /sec . The sea-level acceleration of gravity 
is g^ = 32.20 ft/sec^. 

Atmospheric Model . The atmospheric model used is the US 
Standard Atmosphere, 1976 (Ref. 20). In this model, the values 
of the density are tabulated at discrete altitudes. For intermediate 
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altitudes, the density is computed by assuming an exponential 
fit for the function p (h) . This is equivalent to assuming that 
the atmosphere behaves isothermally between any two contiguous 
altitudes tabulated in Ref. 20. 
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4 . Performance Indexes 

Subject to the previous constraints, different optimization 
problems can be formulated, depending on the performance index 
being considered. The resulting optimal control problems are 
either of the Bolza type or of the Chebyshev type. 

Problem (PI ) . Minimum Fuel Weight. It is required to minimize 
the weight of fuel consumed. Here, the performance index is given 
by 


I = W Q - W 2 (EMI, EM2 ) , 


(20a) 


or 


I = w - w 
0 3 


;em3) 


(20b) 


Problem (P2 ) . Minimum Peak Dynamic Pressure. It is required 
to minimize the peak value of the dynamic pressure. Here, the 
performance index is given by 

I = max (q) , (21a) 

t 

where the dynamic pressure is given by [see (14a)] 


q = (1/2) p v . 


(21b) 


This problem can be reformulated as that of minimizing the integral 
performance index 


J 


q dt 


(EMI, EM2 ) , 


(21c) 


or 
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J = q n dt (EM3 ) , ( 21d) 

J 0 

with n = 8 . 

Problem (P3 ) . Minimum Peak Heating Rate. It is required to 
minimize the peak value of the heating rate at a particular point 
of the aerospace plane, for instance, the stagnation point. Here, 
the performance index is given by 

I = max (Q) , (22a) 

t 

where the heating rate is given by [see (14b) ] 

Q = C/(p/P*) (V/V*) 3 ' 07 . (22b) 


This problem can be reformulated as that of minimizing the integral 
performance index 


J 


f 2 

f Q at 


J 0 


(EMI, EM2 ) , 


(22c) 


or 


J 


(3 

Q n dt 

J 0 


(EM3 ) , 


(22d) 


with n = 8. 

Problem (P4 ) . Minimum Peak Tangential Acceleration. It is 
required to minimize the peak value of the tangential acceleration. 
Here, the performance index is given by 


I = max (a T ) , 


(23a) 
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where the tangential acceleration is given by [see (14c) ] 

a T ~ (Tg e /W)cos(a + 5) - Dg e /W - gsiny. 

This problem can be reformulated as that of minimizing the 
integral performance index 

r 2 


J 


n . , 
a dt 

> 0 


(EMI , EM2) , 


or 

( 3 


J = a£dt (EM3) , 
1 0 


with n = 8. 


(23b) 


(23c) 


(23d) 
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5 . Problem Classification 

In the analyses which follow, the initial path inclination 
Yq is either free or given (Yq ” 0.0 deg); the dynamic pressure 
q is either free or bounded (q £ q ) ; and the tangential acceleration 
a T is either free or bounded (a T £ a ) * Depending on the 
combination of constraints being considered, we have the following 
types of problems: 

(A) Yq — free, q = free, a^ = free, (24a) 

(B) Y n “ free, q < q , a m < a m , (24b) 

(C) Yq = 0*0 deg, q - free, a T = free, (24c) 

(D) Yq = 0.0 deg, q < q^, a T < a^, (24d) 

with 

q u = 1500 lbf/ft 2 , a Tu = 3g e - (25) 

Note that a peak heating rate bound, 

Q £ Q u , ( 26 a) 

with 

Q u = 150 BTU/ft 2 sec, (26b) 

is not imposed because it can be satisfied or nearly satisfied 
indirectly if the dynamic pressure bound is satisfied. 

The ensuing terminology is self-explanatory for the minimum 


fuel weight problem: 
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Problem (P1A) is Problem (Pi), s.t. conditions (A), 

Problem (P1B) is Problem (PI), s.t. conditions (B) , 

Problem (PIC) is Problem (Pi), s.t. conditions (C) , 

Problem (P1D) is Problem (Pi), s.t. conditions (D) . 

A similar terminology is employed for the problems of minimizing 
the peak dynamic pressure, the peak heating rate, and the peak 
tangential acceleration. 

To sum up, for a given aerodynamic configuration and a given 
engine model, there are 16 optimization problems to be solved. 
Since there are three engine models, this leads to a total of 48 
optimization problems to be solved. 
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6. Numerical Results 

Numerical solutions for the optimization problems formulated 
in Sections 4-5 were obtained by means of the sequential gradient- 
restoration algorithm (SGRA, Refs. 21-22) . We recall that SGRA 
is a first-order algorithm which generates a sequence of feasible 
solutions, each characterized by a lower value of the performance 
index being considered. 

While SGRA is available in both primal form (PSGRA) and dual 
form (DSGRA) , the primal form is better suited for hypervelocity 
flight problems; hence, it is employed here. A cross section of 
the solutions obtained is presented in Tables 1-3. For more 
details on the solutions, see Ref. 23. 

6.1. Unconstrained Solutions . First, constraints of Type (A) 
were considered [see (24a)], meaning that Yq is free, the dynamic 
pressure q is unconstrained, and the tangential acceleration a^ 
is unconstrained. For constraints of Type (A) and for engine model 
EMI, SGRA was employed to minimize each of the performance indexes 
of Section 4. Summary results are shown in Table 1, which lists 
the values of the weight of fuel consumed, the peak dynamic 
pressure, the peak heating rate, and the peak tangential acceleration 
for Problems (PlA) , (P2A) , (P3A) , (P4A) . Table 1 also lists the 

initial path inclination y Q , the time duration of each segment 

of the trajectory (t^, T 2^ ' an< ^ = ®2* 

In interpreting the results of Table 1, we stress the 
following concept: even though these results have been obtained 
disregarding the presence of constraints on Yq* *3/ Q* a T r 
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results are acceptable from an engineering point of view only 

if y q , max (q) , max(Q), max(a T ) are in a suitable range. Specifically, 

Yq should be relatively small, otherwise an intolerable burden 

is imposed on the turbojet segment of the trajectory? max(q) 

should be less than 1500 lb/ft 2 ; max(Q) should be less than 

150 BTU/f t 2 sec ; and max(a T ) should be less than 3g e . 

Inspection of Table 1 shows that solution (P1A) is 
not acceptable, mainly because the values of Yq an d max(a T ) 
are excessive? solution (P2A) is not acceptable because the values 
of Yq an< ^ max (a^,) are excessive? solution (P3A) is not acceptable 
because the values of Yq and max(a,p) are excessive? and solution 
(P4A) is not acceptable because the values of Yq/ max(q), max(Q) 
are excessive. In addition, for solution (P4A) , the weight of 
fuel consumed is too large with respect to that of solution (P1A) . 

6.2. Constrained Solutions . Next, Problem (Pi) is considered 
[see (20)], meaning that the weight of fuel consumed is being 
minimized. For Problem (Pi) and for engine model EMI, SGRA was 
employed to obtain minimum fuel solutions for each of the constraint 
combinations (24). Summary results are shown in Table 2, which 
lists the values of the weight of fuel consumed, the peak dynamic 
pressure, the peak heating rate, and the peak tangential acceleration 
for Problems (P1A) , (PlB) , (PIC) , (P1D) . Table 2 also lists the 

initial path inclination Y 0 r the time duration of each segment 
of the trajectory (t^, t 2 ) , and the final time 6 f = 6 2 - 

Inspection of Table 2 shows that solution (PlA) is not 
acceptable mainly because the values of Yq a-nd max(a l j,) 
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are excessive; solution (P1B) is not acceptable because the 
value of Yq is excessive; solution (PIC) is not acceptable 
because the values of max(q), max(Q), max(a T ) are excessive. The 
only acceptable solution is solution (P1D) , albeit with a 1.7% 
increase in fuel weight with respect to that of solution (P1A) . 

6.3. Effect of the Engine Model . Thus far, the only 
acceptable solution is (P1D) , which is obtained by minimizing 
the weight of fuel consumed [Problem (Pi)) in conjunction with 
constraints of Type (D) . Therefore, we consider now Problem (P1D) 
for different engine models (EMI, EM2 , EM3) . 

We recall that engine model EMI is a ramjet/scramjet combination 
with scramjet specific impulse tending to a nearly-constant value 
at large Mach numbers; engine model EM2 is a ramjet/scramjet 
combination with scramjet specific impulse decreasing monotonically 
at large Mach numbers; engine model EM3 is a ramjet/scramjet/rocket 
combination in which the scramjet operates only at M £ 15; at 
higher Mach numbers, the scramjet is shut off and the aerospace 
plane is driven only by the rocket engines. 

Summary results are shown in Table 3, which lists the values 
of the weight of fuel consumed, the peak dynamic pressure, the 
peak heating rate, and the peak tangential acceleration for 
Problem (PlD) and engine models EMI, EM2 , EM3 . Table 3 also 
lists the initial path inclination y Q , the time duration of each 
segment of the trajectory [(t^, t 2 ) for engine models EMI, EM2 
and (t^, t 2 , x^) for engine model EM3] , and the final time [0^ = 0 2 
for engine models EMI, EM2 and 0 f = 0 3 for engine model EM3] . 
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Inspection of Table 3 shows that the constraints on max(q), 
max(Q), max(a T ) are now satisfied or nearly satisfied for 
all engine models. Note that Yq = 0.0 deg for all engine models. 

In percentage of the initial weight, the weight of fuel consumed 
is 34.3% for engine model EMI, 44.3% for engine model EM2, and 
60.7% for engine model EM3 . 

Note that the above results exclude the turbojet segment 
of the trajectory. If one assumes that the weight of fuel consumed 
in the turbojet phase is 5% of the take-off weight, one concludes 
that, in percentage of the take-off weight, the weight of fuel 
consumed is 37.6% for engine model EMI, 47.1% for engine model EM2 , 

and 62.7% for engine model EM3 . 

6.4. Remark . For Problem (PlD) solved in conjunction with 
engine model EM3 , one additional simplification was used. The 
rocket portion of the trajectory, corresponding to M > 15, was 
optimized with the following provision: while the angle of attack 
a(t) is treated as a control, the power setting 3(t) is treated 
as a parameter, which means that $(t) = 0 for M 15. 
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7 . Conclusions 

This paper is concerned, with optimizing the trajectories of 
an aerospace plane. A single-stage-to-orbit (SSTO) configuration 
is considered, and the transition from low supersonic speeds to 
orbital speeds is studied under the following assumptions: the 
turbojet portion of the trajectory has been completed; the 
aerospace plane is controlled via the angle of attack a(t) and 
the power setting £(t) ; the aerodynamic model is the generic 
hypersonic aerodynamics model example (GHAME) . Concerning the 
engine model, three options are considered: (EMI) this is a 

ram jet/ scram jet combination in which the scram jet specific 
impulse tends to a nearly-constant value at large Mach numbers; 
(EM2) this is a ramjet/scramjet combination in which the scramjet 
specific impulse decreases monotonically at large Mach numbers; 
(EM3) this is a ramjet/ scram jet/rocket combination in which, 
owing to stagnation temperature limitations, the scramjet 
operates only at M <_ 15; at higher Mach numbers, the scramjet is 
shut off and the aerospace plane is driven only by the rocket 
engines . 

Under the above assumptions, four optimization problems are 
solved using the sequential gradient— restoration algorithm for 
optimal control problems: (Pi) minimization of the weight of fuel 

consumed; (P2) minimization of the peak dynamic pressure; (P3) 
minimization of the peak heating rate; and (P4) minimization of 
the peak tangential acceleration. 
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The above optimization studies are carried out for different 
combinations of constraints, specifically: initial path 
inclination either free or given (yq = 0) ; dynamic pressure 
either free or bounded (q < 1500 lbf/ft 2 ) ; tangential 
acceleration either free or bounded (a«j» 3g e ) . A peak heating 
rate bound (Q £ 150 BTU/ft 2 sec) is not imposed because it can be nearly 
satisfied indirectly if the dynamic pressure bound is satisfied. 

The effect of the performance index, the constraint type, 
and the engine model on the solutions is studied. From an 
engineering point of view, the most useful solutions are those 
which minimize the fuel weight, while satisfying the 
constraints Y 0 = °, q < 1500 lbf/ft 2 , Q < 150 BTU/ft 2 sec, a T < 3g e . 

The preliminary conclusions are as follows: 

(a) For an aerospace plane governed by GHAME + EMI, the 
SSTO mission requires a weight of fuel consumed equal to 34.3% of 
the weight at the beginning of the ramjet phase, equivalent to 
37.6% of the take-off weight if one includes the turbojet phase. 

(b) For an aerospace plane governed by GHAME + EM2 , the 
SSTO mission requires a weight of fuel consumed equal to 44.3% of 
the weight at the beginning of the ramjet phase, equivalent to 
47.1% of the take-off weight if one includes the turbojet phase. 

(c) For an aerospace plane governed by GHAME + EM3 , the 
SSTO mission requires a weight of fuel consumed equal to 60.7% of 
the weight at the beginning of the ramjet phase, equivalent to 
62.7% of the take-off weight if one includes the turbojet phase. 
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(cl) if one assumes that engine model EM2 is the one closer 
to reality, then the SSTO mission appears to be feasible. 
Obviously, its ability to deliver payloads can be improved via 
progress in the areas of aerodynamic properties and specific 
impulse properties. 

(e) If one assumes that engine model EM3 is the one closer 
to reality, then the SSTO mission appears to be marginal, unless 
substantial progress is achieved in the areas of aerodynamic 
properties and specific impulse properties. Under this scenario, 
alternative consideration should be given to studying the 
feasibility of a two— stage— to— orbit mission (TSTO mission) . 
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Table 1. Unconstrained solutions , : engine model EMI, 

various performance indexes, constraints of Type (A) . 


Quantity 


Problem 



Units 


(P1A) 

(P2A) 

(P3A) 

(P4A) 


(w 0 -w f )/w 0 

0.337 

0.347 

0.357 

0.550 

- 

max (q) 

1540 

999 

1157 

3751 

lbf/ft 2 

max (Q) 

165 

161 

98 

495 

BTU/ft 2 sec 

max(a T )/g e 

9.1 

5.2 

4.0 

1.1 


Y o 

42.0 

50.0 

40.4 

38.3 

deg 

T 1 

34 

54 

48 

144 

sec 

T 2 

409 

475 

731 

704 

sec 

6 f 

443 

529 

779 

848 

sec 

Wf - W 2 and 

6 f = e 2 for 

engine model 

EMI. 
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Table 2. Constrained solutions, engine model EMI, 

minimum fuel weight, various constraint combinations. 


Quantity 


Problem 



Units 


(P1A) 

(P1B) 

(PIC) 

(P1D) 


<w 0 -w f )/w 0 

0.337 

0.340 

0.339 

0.343 

o 

max (q) 

1540 

1112 

1765 

1500 

lbf/ft 

max (Q) 

165 

148 

200 

153 

BTU/ft 2 sec 

max (a T )/g e 

9.1 

3.0 

13.7 

3.0 



42.0 

39.4 

0.0 

0.0 

deg 

T i 

34 

55 

34 

55 

sec 

T 2 

409 

498 

335 

487 

sec 


443 

553 

369 

542 

sec 

Wf = W 2 and 

6 f = e 2 for 

engine model 

EMI. 
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Table 3. Effect of the engine model, Problem (P1D) , 

minimum fuel weight, constraints of Type (D) . 


Quantity 


Engine model 


“Units 


EMI 

EM2 

EM3 


<w 0 -w f )/w 0 

0.343 

0.443 

0.607 

- 

max (q) 

1500 

1425 

1500 

lbf/ft 2 

max (Q) 

153 

157 

110 

2 

BTU/ft sec 

max(a T )/g e 

3.0 

3.0 

3.0 

' 

^0 

o 

• 

o 

o 

• 

o 

0.0 

deg 

x. 

55 

44 

57 

sec 

1 





T 2 

487 

472 

97 

sec 

T 3 

- 

- 

277 

sec 

e f 

542 

517 

431 

sec 

Wf = W 2 and 

6f = 0 2 for 

engine models EMI, 

EM 2 . 


W f = and 

9 f = e 3 f ° r 

engine model EM3 . 
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FIG. 2C • FUEL RATE ( LBF/SEC ) . BETA=1. RAMJET ENGINE. 
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FIG. 3fl . THRUST ( MLBF ) . BETR=1. SCRRMJET ENGINE. 
ENGINE MODEL EMI- 




FIG. 3B • SPECIFIC IMPULSE (SEC). BETfl=l . SCRRMJET ENGINE. 
ENGINE MODEL EMI. 
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FIG. 4R. THRUST ( MLBF ) , BETR=1. SCRRMJET ENGINE, 
ENGINE MODEL EM2 . 
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